Short time decay of the Loschmidt echo 
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The Loschmidt echo measures the sensitivity to perturbations of quantum evolutions. We study 
its short time decay in classically chaotic systems. Using perturbation theory and throwing out 
all correlation imposed by the initial state and the perturbation, we show that the characteristic 
time of this regime is well described by the inverse of the width of the local density of states. This 
result is illustrated and discussed in a numerical study in a 2-dimensional chaotic billiard system 
perturbed by various contour deformations and using different types of initial conditions. Moreover, 
the influence to the short time decay of sub-Planck structures developed by time evolution is also 
investigated. 



I. INTRODUCTION 

Quantum irreversibility studies have become a very ac- 
tive research topic due to a direct connection with quan- 
tum computers and mesoscopics physics |l|, |2| ■ The natu- 
ral quantity for these investigations has been introduced 
by Peres in his seminal paper of 1984 Called later 
as Loschmidt echo (LE) or fidelity, it measures the abil- 
ity of a system to return to an initial state \<f>) after a 
forward evolution with a Hamiltonian H followed by an 
imperfect reversal evolution with a perturbed Hamilto- 
nian H = Hq + SxH' (5x parameterize the strength of 
the perturbation). Thus, it is given by 



M(t) = | (0| exppift] exp[-Lff t] \4>) 



(1) 



(throughout the paper h is set equal to 1). The LE com- 
pares the evolution of an initial state with slightly differ- 
ent Hamiltonians and can distinguish regular and chaotic 
classical dynamics 0, S • 

The LE was recently studied in various chaotic systems 
using several approaches 0, IE S H, 0, E3, E3, 

How- 
ever, few types of decay were discussed in the literature. 
For a very short time, it is straightforward to show that 
the LE has a parabolic behavior M(t) = l-Sx 2 (AH') 2 t 2 , 
with {AH') 2 = [(<j>\H' 2 \<f>) - (<j>\H'\<j)) 2 }. This decay is 
better resembled by the Gaussian function exp[— (t/r) 2 ], 
with characteristic time r = l/(AH'6x). Though this 
regime has experimental relevance |13|. it has not been 
extensively taken into account. 

After this short time decay, a crossover to a per- 
turbation dependent regime was predicted and numer- 
ically observed H H H G3 For very small Sx, in 
which a typical matrix element U of the perturbation 
is smaller than the mean level spacing A, the decay is 
always Gaussian until M(t) reaches its asymptotic val- 
ues M(t — > oo) = Mqo. If U > A, this regime has an 
exponential decay exp(— Ti), with T the width of the lo- 
cal density of states (LDOS). This is usually called Fermi 
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golden rule regime (FGR). When T > A, with A the mean 
Lyapunov exponent of the classical system, a perturba- 
tion independent regime is observed. In this case, the 
decay rate is given by A. Finally, if T exceeds the band- 
width of the perturbation, the LE has a Gaussian decay. 

The properties of the initial state play an important 
role in the behavior of the LE |(J This point can 
be relevant to observe the mentioned regimes. For exam- 
ple, the Lyapunov regime is not displayed if the initial 
state is an eigenfunction of the unperturbed/perturbed 
Hamiltonian [lj. Localized wave packets are needed to 
observe this regime. On the other hand, Zurek has re- 
cently stated that dynamical evolution causes that these 
states develop a sub-Planck structures in phase space, 
and he predicts that these structures enhance their sen- 
sitivity to perturbations |14|. 

In this article we are mainly interested in the short 
time decay of the LE. Disregarding system specific fea- 
tures and the correlations imposed by the characteristics 
of initial state, we show via perturbation theory that r _1 
is given by the width T of the LDOS. In order to see the 
validity of this result in a realistic model, we study the 
characteristic time r in a paradigmatic model of quan- 
tum chaos, a 2-D chaotic billiard perturbed by a contour 
deformation. We regard the influences of the charac- 
teristics of the perturbation and the initial condition to 
this important regime. Therefore, we consider different 
perturbations and various types of initial conditions in 
our numerical study. We find that some perturbations, 
which we call generic, destroy correlations imposed by 
the initial condition giving that t -1 = T. Nevetheless, 
other perturbations do not act in that way and this fact 
produces a slower decay with r _1 < T. In this context, 
we discuss the influence of an initial time evolution of 
the wave packet and the corresponding developed struc- 
tures in phase space in the short time decay of the LE. 
In fact, it is shown that an initial dynamical evolution 
helps to erase the mentioned correlations, with the effect 
of increasing r _1 . If the initial evolved time is smaller 
than the Ehrenfest time, the enhanced decay is described 
entirely by the classical streaching arround the unstable 
manifold given by the Lyapunov exponent [T^. But for 
greater evolved times, the quantum interference lead a 



developed sub-planck structure in phase space and this 
yield that the decay continuos growing. That is, the sen- 
sitivity to perturbations is also enhanced in this case (as 
stated in Ref. 14] ). However , if the perturbation pro- 
duces a decay with r" 1 = T the developed structures in 
phase space do not influence the short time decay. 

The paper has the following structure. Sec. [H] is 
devoted to describe the model system and the various 
shape deformations that we have considered. The paper 
is self-contained with the inclusion of the shape param- 
eter Hamiltonian expansion for a 2-D quantum billiards 
developed in Ref. 16]. In Sec. IIIII the characteristic 
time r is related to the width of the LDOS using per- 
turbation theory. Then, in Sec. IIVI we present the nu- 
merical results. We study the LE for several types of 
initial condition. Starting with the simplest case when 
the initial state is an eigenfunction of Ho, we follow with 
Gaussian wave packets. Finally the initial conditions are 
the evolved Gaussian wave packets in order to study the 
prediction of Ref. 14]. In Sec. |V|we make some final 
remarks. 



II. MODEL SYSTEM: DEFORMED STADIUM 
BILLIARDS 

We use the desymmetrizcd Bunimovich stadium bil- 
liard as a model system to explore the behavior of the LE 
fiof. This paradigmatic system is fully chaotic and has 
great theoretical and experimental relevance 0, 0, 0] . 
It consists of a free particle inside a 2-dimensional pla- 
nar region whose boundary C is shown in Fig. ^ with 
dashed lines. The radius r is taken equal to unity and 
the enclosed area is 1 + tt/4. 

The system is perturbed by boundary deformations 
which preserves the area of the billiard. Deformations 
with different characteristics are chosen in order to un- 
derstand their influences in the short time decay of the 
LE. Fig. QJshows the shape deformations that have been 
considered. The changes of the boundary are parameter- 
ized by 

r(s, Sx) = ro(s) + z(s, Sx) n , (2) 

with s along C, i"o(s) the parametric equation for C, 
and n the outward normal unit vector to C at vq(s) 
(see Fig. IHa)). Case (a), shown in Fig. ^ is 
well described in Ref. 16]. For deformation (b)-(d), 
z(s,Sx) = aSxcos(nNs/P) with a = 0.42, P = 1 + tt/2 
and N=3, 5 and 10. a is chosen in order that the width 
of the LDOS of all deformations are aproximately equal 
in the studied range of Sx . As we shall see in the next 
section, the width of the LDOS is a measure of the mag- 
nitude of a perturbation. 

To solve Eq. (1), eigenvalues and eigenfunctions of 
perturbed and unperturbed system are needed. They 
are obtained using a Hamiltonian expansion of deformed 
billiards which has been recently developed ^(|. F° r 
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FIG. 1: Schematic figure of the system and the various shape 
deformations. Desymmetrised stadium billiard is plotted with 
dashed lines. On continuous lines, deformations of the billiard 
are shown. The curvilinear coordinates used to describe the 
deformations [Eq. 2] is also included. 

Sx < l/ko (fco the mean wave number of the region un- 
der study), the eigenenergies and eigenfunctions of the 
deformed billiard, are connected to the ones of the sta- 
dium by the linearized Hamiltonian expressed in the basis 
of eigenstates at Sx = (from now on we will call <j)^ to 
these eigenstates and the respective eigenenergies) , 

H^{5x)~H^(0) + 5xH'^, (3) 

with 11^(0) = EfiS^v and 

The eigenfunctions and eigenenergies at Sx = are 
obtained using the scaling method 20]. The integral 
above could be viewed as an inner product among the 
wave functions ^jtjf- defined over C. This relation de- 
fines an effective Hilbcrt space in a window Afc « 
Perimeter/ Area The cut-off function Cf M „ = 

exp [-2 (fc 2 - fc2) 2 /(fc Afc) 2 ] restricts the effect of the 
perturbation to states in this energy shell of width fcoAfc. 
It allows to deal with a basis of finite dimension with 
wave numbers around the mean value ko, restricting to 
a particular region Ak of interest. We are considering 
ko = 100 and 2m = 1 (m the mass of the particle) in all 
the numerical calculations presented above. 



III. SHORT TIME DECAY AND THE LOCAL 
DENSITY OF STATES 

Our aim is to characterize the short time decay of the 
LE. As mentioned in the introduction, a simple calcu- 
lation led the short time decay to depend on the initial 
state and on the perturbation . We want to relate the 
characteristic time r of the short time decay with some 
general properties of the perturbation. 

The influence of a perturbation over a quantum system 
could be described by the local density of states (LDOS). 
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The LDOS of an unperturbed eigenstate is denned as 



E 



(M5x)\^)\ 2 5[E - (E v {6x) - 



(4) 

with E v (8x) and <j> v (8x) the energy and eigenfunction 
of the perturbed Hamiltonian [Eq. JSJ]. This function 
shows how the unperturbed states are coupled to the per- 
turbed ones. Because we are not interested in a particular 
state, an implicit average over the unperturbed state /x 
is considered from now on. We have chosen the width 
(dispersion) 

T(Sx) = PvL {E v (Sx) -Ep, Sx)[E v (Sx) -£ M ) 2 , 

(5) 

as a practical measure of this distribution . 

The LDOS exhibits various regimes as a function of the 
strength Sx [2lL l24j . As we shall see, the perturbative 
regime is relevant for our study. Perturbation theory 
(PT) gives the following first order expression for the 
LDOS 



p PT (E,Sx) = S(E) + 



[E v (6x) - EJ\ 



2 S[E-(E v (5x)-E^)]. 

(6) 

Using the definition of the width [Eq. JS}] and Eq. © 
it is straightforward to show 



PT 



This expression works very well for all the perturbations 
of Fig^with strength Sx < 1/ko- 

With these ingredients in mind, let us consider the 
short time decay of the LE. As pointed out previously, it 




is given by M(t) — exp[— (t/r) 2 ] with • 
{c/}\H'\<j)) 2 )Sx 2 . Let be |0) = £ a p 



= ((<f>\H' 2 \ 
the initial state , so 



,H' 2 ]Sx 2 = T 2 T (Sx). 



(9) 



This result is also valid if a. 



assuming that we 



are averaging over several initial conditions. Note that we 
have obtained Eq. JJjJ throwing out all the correlations 
imposed by the perturbation and the wave amplitudes. 
In the next section we will see in what manner these 
correlations influence the decay. 
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FIG. 2: Mean value of the matrix elements H^ v as a function 
of r = (X — v for the deformations (a) (dotted line) and (d) (full 
line). The peak for the deformation (a) is due to the bouncing 
ball orbits. Insets: Image of a piece of the perturbation matrix 
\H'\ shown as a density plot. Fhe left plot corresponds to 
the perturbation label as (a) and the right plot corresponds 
to perturbation label as (d). Clearly, there is a structure in 
the perturbation matrix for pertubation (a) when compared 
to perturbation (d). 



t ~ 2 = E a » H U H >» - (E <v#>*) W (7) 

The perturbation matrix H' v is a banded matrix due 
to the cut-off function Cf^j, (see the insets of Fig. [2J|. 
This band structure is quite generic in realistic systems 
due to the finite range interaction of unperturbed states 
[22l |23|. Inside the band, the matrix elements H' are 
highly fluctuating numbers. At first sight if we ignore 
the system specific features, we can do the diagonal ap- 
proximation of Eq. [7| resulting 

T- 2 = lY / K\ 2 H> 2 „]6x 2 . (8) 

This approximation is also valid for the case in which 
the complex numbers a M behave randomly. Finally, if 
we consider that the \a^\ 2 are nearly constant inside the 
band, 



Each perturbation is characterized by the structure 
and correlations of the matrix elements H',. As an ex- 
ample, Fig. [2] shows the behavior of the mean value of 
H l2 v as a function of the distance r = v — fi for pertur- 
bations (a) and (d) . This function is usually called band 
profile. The cases (b) and (c) are not plotted because 
their band profile are qualitative equal to case (d) . Note 
that perturbation (a) is clearly non-generic due to the 
two important peaks at \r\ ~ 25. This non-universality 
is introduced by the fact that perturbation (a) does not 
connect the bouncing ball states with generic states. In 
the insets of Fig. |3 a density plot of these matrix are 
shown for deformations (a) and (d). We will see in the 
next section that there are correlations between the ma- 
trix elements of the perturbation which are not exposed 
in the band profile but have an important influence to 
the short time decay. 
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IV. NUMERICAL RESULTS 

In this section we study numerically the behavior of the 
short time decay of the LE in the Bunimovich stadium 
billiard perturbed by the contour deformation presented 
in Sec. ^ We consider different types of initial condi- 
tions: Eigenfuction of Hq, Gaussian wave packets and 
evolved Gaussian wave packets . We would like to see 
the range of validity of Eq. [5] for our particular system. 

A. Eigenfuctions of Ho 

The simplest case of the LE is when the initial state 
of Eq. \l\is an eigenstate of Ho. In this case, the LE is 
directly related with the Fourier transform (FT) of the 
LDOS. Then, M(t) is the so called survival probability 
[3 13, defined as 

P{t) = | (^\exp[iHt\ |^) | 2 = \FT[p„(E,Sx)]\ 2 (10) 

As we saw in the previous section, we expect that r for 
P(t) is well described by Eq. P(t) is the mean value 
of the survival probability over several initial states. We 
have numerically observed that this is the case for the 
perturbations (b-d). These results are shown in Fig. [3] 
Note that t" 1 for perturbations (c) and (d) are not plot- 
ted because the results are the same as (b). The widths 
T(6x) and Tpt(Sx) are equal for all the perturbations. 
For perturbation (a), t _1 rs 0.85r(5a;) which is directly 
related to the structure of the perturbation matrix (see 
Fig. |2J imposed by the bouncing ball states El • Similar 
attenuation was observed in Ref. ^} f° r perturbation 
(a) in the FGR regime. In this case, the decay rate is 
given by 0.5T(Sx) insted of the expected value T(8x). 

Fig. ^summarizes the behavior of P(t) in the desym- 
metrized stadium billiard. We have taken an average 
over 100 initial states. In this figure, the results for per- 
turbation (d) of Fig. 1 are shown. Other perturbations 
led to the same qualitative results. For a small pertur- 
bation strength (top curve of Fig. [2J we observe the 
Gaussian short time decay and after that an exponential 
decay with a decay rate given by the width of the LDOS 
0, llOj . For large perturbation strength (bottom curve 
of Fig. 0} the exponential decay with decay rate given 
by T(Sx) is not observed . This is due to the constraint 
imposed by the asymptotic value M^,. In this case the 
decay is completely Gaussian. An important point is that 
the asymtotic « A/T(Sx). 

B. Localized Gaussian wave packets 

The decay of the LE for localized Gaussian wave pack- 
ets has been widely studied in the literature. Most of 
the previous works consider this case 0, El IS El El- 
The predicted crossover from a perturbation dependent 
regime to the Lyapunov regime has been shown for these 
classically adapted initial conditions. 




5x 



FIG. 3: t -1 of the survival probability as a function of 5x 
when the system is perturbed by deformations (a) and (b). 
The width of the LDOS T(5x) and its perturbative evaluation 
Tpt(Sx) calculated using Eq. (6) are also plotted. 
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FIG. 4: P(t) for the desymmetrized stadium billiard per- 
turbed by deformation (d) of Fig. 1. The values of the 
strength of the perturbation is Sx = 0.015 for the top curve 
and 0.05 for the bottom curve. The dashed lines correspond 
to the Gaussian decay exp[— (F(Sx)t) 2 } and the exponential 
decay exp[— (F(Sx)t)] is plotted with a dotted line. 

We discuss here the short time decay of the LE for 
this particular type initial conditions in the stadium bil- 
liard perturbed by deformations presented in Fig. 1 . Wc 
compute M(t) for initial Gaussian wave packets, 

0(r) = {na 2 )^ 2 exp[ip .(r-r Q )-i\r-r \ 2 /a 2 } (11) 

with \p \ = k = 100 and a = 0.16. An average over 50 
initial states was taken. The direction of the momentum 
Po/|pb| an d the center of the wave packet r*o are cho- 
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sen randomly. As expected, the short-time decay is well 
described by a Gaussian function exp[— (t/r) 2 ]. Fig. |S] 
shows the behavior of r _1 as a function of the strength 
Sx for all the perturbations under study. This type of 
initial conditions imposes correlations so that the sec- 
ond sum of right hand side of Eq. Q does not vanish. 
The non-universalities of each perturbation are clearly 
exposed in r. Note that these differences are not seen 
in the width of the LDOS nor in band profile of the ma- 
trix H' . We have considered perturbations with greater 
number of oscillations of the boundary (N > 10) and we 
find r = l/T(Sx) for all of them (26|. 
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FIG. 5: r -1 of M(t) as a function of Sx when the system is 
perturbed by deformations (a), (b), (c) and (d). The initial 
states are the Gaussian wave packets of Eq. (11). T(Sx) is 
also plotted in solid line. 



C. Evolved Gaussian wave packets 



has been evolved a time To is given by 

r- 1 (T )^r- 1 (0)exp( A(r ° 2 ' rc) ), (12) 

with t(0) the characteristic time of the short time de- 
cay for initial states that have not been evolved, A the 
Lyapunov exponent and T c the mean time for the first 
collision with the boundary. 

We consider the influence of the developed structures 
in phase space in the short time decay of the LE. So, M{t) 
is computed for the same initial conditions [Eq. (flip ] of 
the previous section but an unperturbed initial evolution 
during a time To is applied. Fig. El shows t _1 (To) as a 
function of the strength Sx for perturbation (a) and with 
T = 0,0.025,0.05, and 1. Note that t _1 (T ) increases 
with larger To, for all perturbation strenghts. This fact 
clearly points out that an initial evolution enhances the 
sensitivity to this particular perturbation. We have ob- 
served that for To > 0.4, r _1 converges to the width of 
the LDOS. Same behavior is shown when the system is 
perturbed by deformations (b) and (c). However, when 
the perturbation destroys the correlations imposed by the 
initial wave packet which implies that r = 1/r, the short 
time decay is not affected by an initial evolution. 
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In a recent paper 14], Zurek shows that dynamical 
evolution of initial Gaussian wave packets in classically 
chaotic systems, produces finer and finer phase space 
structures which saturates after the Ehrcnfcst time with 
a sub-Planck scale. More important, he predict that this 
sub-Planck structures enhances the sensitivity of a quan- 
tum state to an external perturbation. A numerical study 
in a time dependent one dimensional model agreed with 
this assertion |27[. but other studies reached opposite 
conclusions |15L |28|. In Ref. ^5| it is showed an en- 
hanced decay of the LE of evolved wave packets but this 
acceleration is fully described by the classical Lyapunov 
exponent and it is not due to the sub-Planck structures. 
More specifically, it is pointed out that the characteristic 
time t(Tq) of the short time decay for a wave packet that 



FIG. 6: t -1 of M(t) as a function of Sx for perturbations (a). 
The initial states are Gaussian wave packets that have been 
evolved a time To. T(Sx) is also plotted. 

In order to see if these increments are fully described 
by Eq. 1121 and due to that has classical nature, in Figd 
it is showed the behavior of r _1 (To)/r _1 (0) for several 
preparation time To for perturbation (a) with strenght 
5x = 0.04. It is clearly observed that Eq. El works 
well for T < 0.025. Note that the Ehrenfest time T E = 
0.025. So, the enhacement of the short time decay is 
fully explained with the classical incstability given by the 
Lyapunov exponent for evolved times T smaller than the 
Ehrenfest time. However, for larger times r _1 (To) is also 
growing. 
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A qualitative picture of acceleration of the short time 
decay for evolved states is the following. Before the 
Ehrenfest time the wave packet is streaching around 
an unstable manifold and just after that times starts 
the quantum interference which lead a sub-Planck struc- 
tures. At that times a small part of the Wigner function 
presents a sub-Planck structure. This region grows with 
time and it seems to be the reason of the accelerating 
decay. Note that for the saturation times the sub-Planck 
structure is all arround the available phase space. 
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FIG. 7: T _1 (Tb)/r _1 (0) as a function of the preparation time 
To for strenght Sx = 0.04. The stadium is perturbed by de- 
formation (a). The prediction of Ref. fTEf (Eq. I12B is ploted 
in solid line. 



V. FINAL REMARKS 

We have studied the short time decay of the LE in 
a 2-D chaotic billiard. The system was perturbed by a 
contour deformation. Different perturbations were con- 
sidered in order to develop the influences of their char- 
acteristics in the behavior of the LE. Moreover, several 
types of initial conditions have been used and how they 
affect the LE have been examined. 

Our findings are the following. For non-localized ini- 
tial states and if the system is perturbed by a generic 
deformation, the characteristic time t of the short time 



Gaussian decay is given by the inverse of the width T 
of the LDOS. If semiclassical features are exposed in the 
matrix elements of the perturbation, we have obtained 
that r _1 < r. For highly localized initial states, cross 
correlation between wave amplitudes are important and 
this is exposed with the fact that r _1 < T. When the 
perturbation destroys such correlations, the characteris- 
tic time t -1 exhibits its maximum value T. 

We have discussed the prediction of Zurek [l4| which 
stated that an initial dynamical evolution of semiclassi- 
cal wave packets lead a sub-Planck structures in phase 
space and this enhances its sensitivity to perturbation. 
We found that in the cases in which the perturbation 
does not destroy the correlation mentioned before an ac- 
celerated decay is observed. As a function of the prepa- 
ration time To, we have observed two regimes. For To 
smaller than the Ehrenfest time, the enhanced decay is 
described entirely by the classical Lyapunov exponent as 
pointed out in Ref. jl5| |. However, for larger To in which 
the quantum interference lead the sub-Planck structures 
the enhancement of decay is also observed. 

A final point is worth commenting. We have shown 
that the short time decay of the LE has a Gaussian be- 
havior and for certain perturbations the characteristic 
time is given by the width of the LDOS. The LDOS of 
some systems 0,0,HEI13 presents a region in which its 
width is independent of the perturbation. We note that 
these results could be of importance for the understand- 
ing of the measure of the LE in recent NMR experiments 
Although that system consists of many interact- 
ing nuclear spins, the results are in accordance with the 
former. That is, due to imperfections in the reversed evo- 
lution the LE in the MNR experiment shown a Gaussian 
attenuation, M(t) = exp[(i/r) 2 ], with t depending on 
the small non-inverted interaction (which corresponds to 
the perturbation SxH' of Eq. , and in a range of small 
perturbation, the characteristic time r do not depend on 
it. 
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Our system do not exhibit such a region due to inclusion 

of the cut-off function. This region can be observed for 

large enough 8x but in this case we can not deal with the 

numerics. 



